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Abstract 

We study a Bose-condensed gas at finite temperature, in which the particles of the 
condensate and of the thermal cloud are constrained to move in a plane under radial 
harmonic confinement and interact via strictly two-dimensional collisions. The cou- 
pling parameters are obtained from a calculation of the many-body T-matrix and 
decreases as temperature increases through a dependence on the chemical potential 
and on the occupancy of excited states. We discuss the consequences on the conden- 
sate fraction and on the density profiles of the condensed and thermal components 
as functions of temperature, within a simplified form of the two-fluid model. 
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1 Introduction 



Recent experiments have realized quasi two-dimensional (2D) cold Bose sys- 
tems by tuning the anisotropy of the trapping potential [1,2] and have stim- 
ulated interest in studying Bose-Einstein condensation in a low-dimensional 
regime. An important difference between the macroscopic 2D fluid and the 3D 
one is that at finite temperature phase fluctuations in 2D destroy the long- 
range order, in agreement with the Mermin- Wagner theorem [3] . A Kosterlitz- 
Thouless transition to a superfluid state still occurs on cooling through the 
binding of vortices of opposite vorticity [4] and leads to an algebraic decay 
of the one-body density matrix. However, in a trapped 2D fluid the modifi- 
cations of the density of states caused by the confining potential allow a true 
condensate to exist even at finite temperature. 
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Another important consequence of lowered dimensionality is that the T-matrix 
for two-body collisions in vacuo at low momenta and energy, which should 
be used to obtain the coUisional coupling parameters to lowest order in the 
particle density, vanishes in the strictly 2D limit as the s-wave scattering length 
becomes larger than the width of the axial trapping [5,6]. It is then necessary 
to evaluate the scattering processes between pairs of Bose particles by taking 
into account the presence of a condensate and a thermal cloud through a 
many-body T-matrix formalism [7~9] . 

This formalism has already been used in a number of studies of low dimensional 
Bose gases, dealing in particular with phase fluctuations and the Kosterlitz- 
Thouless transition in a variational approach [7,8,10,11], with a mean-field 
evaluation of the breathing-mode frequency in a trapped ID gas [12], and 
with the equilibrium density profile of a 2D condensate [9,13]. In the limit 
of zero temperature the condensate-condensate coupling parameter has been 
related to the two-body T matrix by considering that, when two particles 
in the condensate collide at zero momentum, they both require an energy 
equal to the chemical potential /i to be excited out of the condensate [9]. 
With increasing temperature the population of the excited states becomes 
non-negligible and a microscopic theory which also takes into account the 
depletion of the condensate is required [14] . An additional coupling parameter 
to describe the scattering processes between an atom in the condensate and 
an atom in the thermal cloud has been introduced by Stoof and co-workers 
[10] through the two-body T-matrix at energy — 

The present paper is the finite-temperature extension of the work by Tanatar et 
al. [13], who studied the cross-over from a 3D to a 2D regime in the equilibrium 
density profiles at zero temperature. Here we limit ourselves to the strictly 2D 
case, where the scattering length has become larger than the axial thickness of 
the cloud. We give a semi-analytical expression for the many-body T-matrix 
elements corresponding to the condensate-condensate and to the condensate- 
thermal cloud coupling parameters {g2 and gi, respectively) and we use them 
in a two-fiuid model [15] to evaluate the condensate fraction and the density 
profiles of the bosonic gas at increasing temperature. Congruently with the 
work of Tanatar et al. [13] we focus on the experimental parameter of Gorlitz 
et al. [1] for the case of largest anisotropy. 

The paper is organized as follows. Section. 2 summarizes the general treatment 

of the many-body T-matrix for two-body scattering processes in the presence 
of a condensate and a thermal cloud. In Sec. 3 we evaluate a simplified form 
of the two-fiuid model for a 2D Bose gas at finite temperature. Finally, we 
present and discuss our results in Sec. 4. 
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2 Coupling pcirameters in a mixture of condensed and thermal 
bosons 



The scattering processes between pairs of atoms in a gas consisting of a Bose- 
Einstein condensate and a thermal cloud are described by the many-body 
T- matrix T^^{E) as a function of the energy E. The coupling parameters 
are given by the matrix elements (k'| T^^^[E) |k) = T*^^(k, k', K; E'), taken 
in the limit of zero energy and momenta. Here k and k' are the incoming and 
outgoing relative momenta of the pair of center-of-mass momentum K. We 
shall consider only the condensate-condensate and condensate-thermal cloud 
couplings and neglect the scattering between thermally excited atoms in the 
present case of a dilute gas. 

Before discussing the many-body T-matrix, however, we shall first recall the 
behavior of the two-body T-matrix that describes collisions between pairs of 
particles in vacuo. We shall then discuss how the two T-matrices are related 
in the appropriate limit for a two-fluid system. 



2.1 The two-body T matrix 



The two-body T-matrix is the solution of the Lippmann-Schwinger equation, 

(k'lr^^(^) |k) = (k'|y(|ri-r2|)|k) 

+ E (k'l ^(|ri - r^l) |q) (q| T'^{E) |k) , (1) 

q E- le^ 

with T^(|ri — r2|) being the interparticle potential. The center-of-mass energy 
for the pair of atoms is E and each atom of mass m has single-particle excita- 
tion energy Cq = /i^g^/2m, as collision takes place in free space. In Eq. (1) the 
collision is described by a single-loop interaction between the two atoms plus 
contributions involving all possible transition routes from state |k) to state 
|k') via intermediate states |q). 

In the case of hard-disk potential of strength Vq = 47r^^/m and in the dilute 
limit ka, k'a -C 1 where a is the 2D scattering length, the solution of Eq. (1) 
is [10] 

In this case the T-matrix is independent of the momenta, but depends on the 
logarithm of the energy and vanishes as E approaches zero. Therefore, the 
presence of the surrounding gas must be taken into account in the coUisional 
processes. Following the proposal of Morgan et al. [16], this is done by setting 
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E — —2// for the mutual scattering of two condensate particles, this being the 
energy required for them to be excited out of the condensate. By a similar 
argument Al Khawaja et al. [10] set E — —fi for the scattering between a 
boson in the condensate and a boson in the thermal cloud. 



2.2 The many-body T-matrix 



According to the above argument, a first approximation for the many-body 
T-matrix at zero momenta and energy is 

rf^(0,0,0;0) = r2^(0,0,0;-n/.) (3) 

with n = 1 for condensate-thermal cloud scattering and n = 2 for condensate- 
condensate scattering. Further many-body effects can enter the scattering pro- 
cesses attended by the presence of a condensate from attributing a Bogoliubov 
spectrum to the intermediate states [7]. At low momenta the many-body T- 
matrix is the solution of the integral equation, 



(k'l T^^{E) |k) = (k'l V{\r^ - r^\) |k) + ^ (k'| V(|ri - r^l) |q) 

q 

1 + 2N{u^) 



X- 



E-2nw^ 



q|T^^(S) |k) 



(4) 



where hw^^ ~ eq -|- are the Bogoliubov excitation energies in the Hatree ap- 
proximation II <^ eq and N{u!q) ~ {exp[/3(eq-|-//)] — are the corresponding 
population factors. For a contact interaction potential Eq. (4) yields 



r^^(o,o,0;£;) 



1 _ ^ 1 -t- 2iV(a;q) 



Vo 



q E-2hw^ 



(5) 



As discussed by Stoof and Bijlsma [7,8], Eq. (5) is still affected by an ultraviolet 
divergence, which can be remedied by replacing Vq in favor of the two-body 
T-matrix. The final result is 



0, 0; 0) = T2>, 0, 0; -n/i) 



<2B( 



l + r2^(0,0,0;-n//)^ 



q 



, (6) 



where 7^-^(0, 0, 0; —njj) is given by Eq. (2) with E = —nji. 



2.3 Calculation of coupling parameters 



The sum in the RHS of Eq. (6) can be evaluated analytically by replacing the 
sum over intermediate states by an integral over momentum. Setting iV(a;q) = 
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E^iexp(-s/3;k^q), we get 



m ~ 



where Ei{x) is the exponent-integral function. In the asymptotic low-temperature 
regime ^ 1) this thermal-population term can be approximated by 

oo 

J2Ei{-sPf^) ^ (/3/^)~'ln[l - exp(-/3/x)]. (8) 

s=l 

A numerical illustration is given in Fig. 1 for values of the system parameters 
appropriate to the experiments on ^^Na by Gorlitz et al. [1] (see Sec. 4). 

In summary, the coupling parameters in the 2D Bose gas are given by 

MB ^ Anhym 

\n\Ahy{nmfia^)\-2j:T=iEii-s(3i2)' 

with n=2 for collisions between pairs of condensate atoms and n=l for colli- 
sions between an atom in the condensate and a thermally excited atom. These 
parameters depend on temperature both through the chemical potential and 
through the excited-state population factor given in Eq. (7) and asymptoti- 
cally approximated by Eq. (8). If this population factor is dropped, one obtains 
from Eq. (3) the corresponding "two-body" coupling parameters as 

= . (10) 

In the next Section we evaluate the chemical potential and hence the coupling 
parameters in Eqs. (9) and (10) through a self-consistent evaluation of the 
density profiles in a two fiuid-model. 



3 Equilibrium properties in the two-fluid model 



Let nc{r) and Undr) be the particle density profiles for the condensate and for 
the thermal cloud in a 2D Bose gas which is radially confined inside a isotropic 
planar trap described by the external potential Vextir) = The two- 

fluid model [15,17] combines a solution of the Gross-Pitaevskii equation for the 
condensate with a Hatree-Fock model of the thermal cloud, which is treated 
as an ideal gas subject to an effective potential K//(^)- In the present case, 

Veff(r)^Vext(r) + 2ginc. (11) 

As already noted, we are neglecting the collisions between pairs of bosons 
belonging to the thermal cloud. It is well known that a full numerical solution 
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of the Gross Pitaevskii equation can be avoided when the kinetic energy term 
in it can be neglected [18]. This yields the condensate density in the Thomas- 
Fermi approximation as 

n,{r) = (l/^2)[/i - Kxt(r) - 2^in„e(r)]^(/i - Kxt(r) - 2g,nnc{r)). (12) 

Since we are interested in examining qualitative behaviors rather than in at- 
taining high numerical accuracy, we have adopted Eq. (12) for the condensate 
density and found that discontinuities occurring in the density profiles at the 
Thomas-Fermi radius can be eliminated by the simple expedient of introduc- 
ing momentum cut-off in the expression for nnc{r). That is, we calculate findr) 
from 



(13) 



where we take po = ^J2ri^f^n^c [19]. The model is then evaluated by solv- 
ing self-consistent ly Eqs. (11)-(13) together with the condition that the areal 
integral of nc(r)-|-n„c(^) should be equal to the total number N of particles. 

Before turning to a presentation of our numerical results, let us point out the 
wide range of validity of the present model. A mean-field treatment is valid 
when the diluteness condition UcO^ <^ 1 holds and if the temperature of gas 
is outside the critical region. With regard to the thermal cloud, no significant 
differences have been found between the predictions of the Hatree-Fock and 
Popov approximations in the regime ficO^ <^ 1, except at very low temperature 
where the thermal cloud is becoming negligible [20]. 



4 Results and discussion 



For a numerical illustration, we have taken the values of particle number, the 
radial trap frequency, and the scattering length as appropriate for ^^Na atoms 
in the experiment of Gorlitz et al. [1] {N ^ 5 x 10^, cu = 188.4 Hz, a = 2.8 
nm). Whereas in their experiment collisions are in a 3D regime, we focus on a 
strictly 2D regime that could be reached experimentally by increasing either 
the trap anisotropy parameter or the scattering length. 

In Fig. 2 we compare the temperature dependence of the coupling parameters 
g^^ (long dashes for n=2 and short dashes for n=l) with that of the coupling 
parameters g^f^ (full and dotted lines, respectively). It is evident that the 
many-body screening of the interactions due to the occupancy of excited states 
is quite large and rapidly increasing with temperature. 

Such many-body screening has, however, very little effect on equillibrium prop- 
erties of the gas for our choice of system parameters. Figure 3 reports the 
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condensate fraction Nq/N for the g^^ model as a function of temperature, in 
comparison with that of an ideal Bose gas at the same values of the system 
parameters. As is well known, the transition temperature and the condensate 
fraction are lowered by the interactions. However, the g^^ model gives results 
that are practically the same as the g!^^ one. 

Finally, the panels in Fig. 4 show the evolution of the density profiles for the 
condensate and for the thermal cloud with increasing temperature from near 
absolute zero to the critical temperature ~ O.QSTq with Tq = {\/6/7r)h(jjVN 
being the critical temperature of the ideal gas [21]. The results in Fig. 4(a) are 
in good agreement with those of Tanatar et al. [13], except that the tails of the 
profile are missed in the Thomas-Fermi approximation. Again the screening of 
collision from the occupancy of the excited states is very small and becomes 
barely visible at T 0.75 %. 

In conclusion, the use of Eq. (3) to describe the many-body effects in two- 
body scattering processes in 2D Bose-condensed gas appears to be very good in 
regard to equillibrium properties at large values of particle number. A decrease 
in the number of particles lowers the chemical potential and may lead to 
observable effects for 10^. 
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Figure captions 



Figure 1. The correction term A{(3fj,)=J2^=iEi{—nPn) from excited-state 
occupancy (dashed hne) and its approximate form from Eq. (8) (sohd hne) as 
functions of P/j,. In the inset a zoom of the region 0.2 < < 1 is shown. 

Figure 2. Interaction strengths (in units of ti? /m) as functions of temperature 
T (in units of ideal-gas critical temperature Tq). 

Figure 3. Condensate fraction Nq/N as a function of temperature T (in 
units of To) from the g^^ model and (dashed line) compared with the non- 
interacting gas (solid line). 

Figure 4. Density profile of the condensate (solid line) and the thermal cloud 
(dashed line) in the g^^^ model (in units of N/a\^ with aho = {h/muY^'^) 
as functions of radial distance r (in units of aho) at various values of the 
temperature (T/Tq^O, 0.3, 0.5, 0.8 and 0.95 from (a) to (e)). 
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0.5 1 1.5 2 2.5 3 

Fig. 1. The correction term A{[3^)=Y^^^Ei{—s[3^) from excited-state occupancy 
(dashed hne) and its approximate form from Eq. (8) (sohd hne) as functions of fiji. 
In the inset a zoom of the region 0.2 < [3ijl < 1.2 is shown. 
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Fig. 4. Density profile of the condensate (solid line) and the thermal cloud (dash( 

line) in the model (in units of N/a^ with aho = (^/mu;)^/^) as functions 
radial distance r (in units of aho) at various values of the temperature (T/To= 
0.3, 0.5, 0.8 and 0.95 from (a) to (e)). 
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